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Abstract The global structure of the current flows in pulsar magnetospheres 
is investigated, with rough calculations of the elements in the magnetospheric 
circuit. It is emphasized that the potential of critical field lines is the same as that 
of interstellar medium, and that the pulsars whose rotation axes and magnetic 
dipole axes are parallel should be positively charged, in order to close the pulsar's 
current flows. The statistical relation between the radio luminosity and pulsar's 
electric charge (or the spindown power) may hint that the millisecond pulsars 
could be low-mass bare strange stars. 
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1 INTRODUCTION 

Since Hewish et al. (1968) discovered the first radio pulsar, more and more magnetospheric 
models for pulsars have been proposed to explain their observed phenomena. The vacuum inner 
gap model, on the one hand, suggested first by Ruderman & Sutherland (1975, hereafter RS75) 
depends on enough binding energy of charged particles on the pulsar surface. The space charge- 
limited fiow model (e.g., Arons & Scharlemenn 1979; Harding & Muslimov 1998), on the other 
hand, is without any binding energy. An outer gap near light cylinder has also been proposed 
(e.g., Cheng, Ho, & Ruderman 1986), the existence of which may also reflect strong binding 
of particles on pulsar surface (Xu 2003a). However, as the observational data in radio, optical. 
X-ray, and 7-ray bands accumulates, there are still a great number of puzzles to be solved 
(e.g., Melrose 2004). Nevertheless, it is obvious that the models depend on the nature of pulsar 
surfaces, and one may conclude on the interior structure of pulsars (e.g., either normal neutron 
stars or bare strange stars, see, e.g., Xu 2003b), which is almost impossible via calculations in 
supranuclear physics, via investigating pulsar emission models. 

The study of pulsar magnetosphere is essential for us to understand various radiative pro- 
cesses, and thus observed emission in different bands. Goldreich and Julian (1969) argued that 
a pulsar must have a magnetosphere with charge-separated plasma and demonstrated that a 
steady current would appear if charges can flow freely along the magnetic field lines from the 
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pulsar surface. Sturrock (1971) pointed that such a steady flow is impossible due to pair for- 
mation, and suggested that a simple electric circuit of the pulsar magnetosphere would be a 
discharge tube connected with an electromotive source. RS75 proposed that sparking process 
takes place in a charge-depletion gap just above the pulsar surface. The sparking points drift 
due to E X B, which can naturally explain drifting subpulse phenomena observed in radio band. 
Based on the assumption that the magnetosphere has a global current loop which starts from 
the star, runs through the outer gap, the wind and the inner gap, and returns to the star, 
Shibata (1991) proposed a circuit including an electromotive source connected in series with 
two accelerators (the inner and outer gaps) and the wind. Providing a fully general relativistic 
description, Kim et al. (2005) studied the pulsar magnetosphere and found that the direction 
of poloidal current in neutron star magnetosphere is the same as that in black hole. 

In this work, assuming that the critical magnetic-field lines are at the same electric poten- 
tial^ as the interstellar medium (ISM, Goldreich & Julian 1969), we propose a circuit model for 
pulsar magnetosphere. The corresponding relation of the elements between magnetosphere and 
circuit is as follows. (1) The total inner radiation of the pulsar corresponds to an electromotive 
power parallel connecting with a capacitor. Note that the resistance of power is negligible due to 
perfect conductivity of the star. (2) The inner gaps which include the inner core gap (ICG) and 
the inner annular gap (lAG) (Qiao et al. 2004a, 2004b) correspond to a parallel connection of a 
resistor and a capacitor. When a spark takes place in gap, it can be represented by a resistor; if 
there is no spark, the voltage on the gap is so high that it can be described by a capacitor. (3) 
The outer gap and the pair-plasma wind correspond to a series connection about a inductor and 
a parallel connection of a resistor and a capacitor. Considering a parallel pulsar whose rotation 
axis is parallel to the magnetic axis, we find that the pulsar should be positively charged on 
the surface for the necessarian of a close circuit. The total electric field along field lines, E^\, 
in the magnetosphere is then composed of two components: that due to charge- departure from 
the Goldreich- Julian density, and that due to the charges of pulsar. 

This paper is arranged as follows. The model is introduced in §2. The total charges of pulsar 
are estimated in §3. In §4, we discuss about the charges of low-mass strange stars and show 
evidence for low-mass millisecond pulsars by observational data. Conclusions are presented in 
§5. 

2 THE MODEL 

As shown in Fig.l, the foot points of line a, b, and rotation axis on the star surface are assigned 

as A, B and P. Point P is also the magnetic pole of star. We assume that the potential of the 
critical field line b equals to that of ISM at infinity, (pB = 0', otherwise, a close electric current 
in the two regions (i.e.. Region I: with boundary lines labelled "a", and Region II: that between 
lines "a" and "b") of open field lines is impossible^. Then for a parallel pulsar, the potential 
(pi < (within Region I) and (pu > (within Region II) in the regime of 0b = 0. Therefore, the 
negatively charged particles should flow out along the open magnetic field lines within Region 
I from star, but positively charged particles flow out through Region II. 

The inner gaps (including ICG and lAG in this paper) and an outer gap may work in a 
magnetosphere^. There is no current in the circuit until a spark takes place in the inner gaps, so 

^ We choose the potential of the interstellar medium to be zero, 0ism = 0, in this paper. 

^ For parallel pulsars, electric current flows outwaxd in the two regions if one sets the potential of the 
polar line to equal to that of ISM, (pp — 0, since all the potentials of open-field lines, except the polar 
line, arc greater than zero. Also one can see that the current flows inward in those two regions if one 
chooses the potential of the last open-field lines to be zero, (pA = 0. Current flows can not be closed in 
both these cases. 

* The inner annular gap and the outer gap might not exist simultaneously. 
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Fig. 1 Sketch of the magnetic field fine, electric field line, and magnetospheric charge 
distribution for a parallel pulsar. Line a is the last open magnetic field line. Line b 
is the critical field line. Region I: with a boundary of lines labelled "b" , and Region 
II: between lines "a" and "b" . If the potential of line 6 equals to that of interstellar 
medium, then the potential within Region II is positive but negative within Region I. 



these gaps could be simulated by parallel connected a resistor and a capacitor, which connect 
with other parts in series (as shown in Fig. 2). 

In Fig. 2, the power of star could be equivalently modelled by voltage V and capacitance 
Cstar- The star is magnetized and possesses an interior electric field, which satisfies 

E+ xB = 0, (1) 

c 

where VL is the angular velocity of star rotating around the dipole rotation axis, which relates 
with rotating period of star P by P = 2n/Q,. The magnetosphere of a rotating isolated-pulsar 
is thus generally concluded to be powered by an electric source with certain potential drop 
between the polar angle 9^ and 9 of 

(j)= (sin^ 9 - sin^ ^b) ?^ 3 x 10^^P|Pi2P~^ (sin^ 9 - sin^ 6>b) volts, (2) 

if the pulsar is assumed to be magnetized homogenously, where Pg = P/(10^ cm), P12 = 
P/(10^^ G). According to the equation of dipolar field line, r = rd sin^ 6* (rj is the maximum 
polar radius), and the polar angle of null surface = cos~^(±l/'\/3), one can obtain sin^ 9b = 
(2/3)3/2p/PL , where the radius of hght cyhnder Pl = c/VL ^ cP/(27r). 
The capacitance of the star with a radius P is, in cgs units. 



Cstar = 9 X lO^^^P farads. 



(3) 
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Fig. 2 An equivalent circuit of pulsar magnetosphere. An electromotive source con- 
nected with other parts in series corresponding to the ICG, lAG, the outer gap and the 
pair-plasma wind. The current is determined by the power voltage V. "i?" stands for 
plasma resistance effects if a spark happens, mainly the losses of relativistic particles. 
"C" represents the gap capacitance effects if there is a potential drop on the gap. "L" 
describes the electromagnetic characteristic of flow. 



The electric flow process in a magnetosphere can be equivalently mimicked by an antenna cable 
with capacitance and inductance. They could be equivalently described by the inductance Lflow 
and Cfiow, the values of which can be estimated as 

Cflow ^2^eoKln-)-' (4) 
n 

and 

iflow ^ — iln— (5) 

ZTT ri 

respectively; where I, ri, and r2 are the length, inner radius, and outer radius of the cable, 
respectively. The passage length of the electric current / ~ y/3/2ri^ — 5.8 x lO^P m. Near the 
light cylinder, the ratio r2/ri ^ (rL -I- rL/V2)/rh = 1.7; but on the stellar surface, the ratio 
r2/ri — (3/2)'^/'* = 1.4. Both them are not dependent on P, B, and other parameters, and we 
choose thus simply r2/ri = 1.5 in this paper. One therefore comes to 

Cflow 8.0 X 10~^P farads, (6) 

and 

iflow ~ 4.7P henries. (7) 

In case enough binding energy of charged particles on the stellar surface, an RS-type (RS75) 
vacuum gap should exist near polar cap, which can be the equivalent of a capacitor of two 
parallel slabs, with 

Crs = neorl/h = 5.8 x IQ-^RlPh^^ farads, (8) 

where rp = _Rsin Op = 1.45 X lO'^Pg/^pi/^ cm is the radius of polar cap, the gap height h is a 
model dependent parameter, /13 — h/{10'^ cm). 
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For the curvature-radiation-induced and the resonant inverse-Compton-scattering-induced 
cascade models, the gap heights h could be (e.g., Zhang, Harding, & Muslimov 2000) 

h„ = 54pl/'B-2^/'p^/' m, (9) 

and 

/lies = 279p^/X2"^'^'/' m, (10) 

respectively, where p is the radius of curvature of field line, pe = p/(10^ cm). 

If the field lines which cross the light cylinder can only not co-rotate, the active region of the 
star is then the polar cap, from polar angle to 6'p ~ (cP) = 1.45 x 10~^(Rq/P)~^/^ . 

In this case, the potential difference of the electrical source is therefore 

1/ap = -6.58 X lO^^iJ^BiaP-^ volts. (11) 

Since the dominate source of rotation energy dissipation is through .Rwind) we can estimate 

-^wind 

R^ind^-^ = nM^^RiBf2Pi5^p-' ohms, (12) 

where the rotation loss rate Erot = — 47r^/_P/P^, I ~ lO'^'^MiRl g-cm^ is the moment of 
inertia for a neutron star with mass ^ MxMq and radius ~ i?6 x 10 km, the period derivative 
A 5 = |-P|/10^^'''. Compared with i^windi the stellar resistance iistar is negligible due to the 
perfect conductivity of the star. 

The potential drop of outer gap, where no spark happens, corresponds to a resistor, which 
is presmned to be combined with the wind dissipation as a total one -Rwind hi Fig. 2. In the mag- 
netodipole radiation model, the filed B and the spindown P is connected by (e.g., Manchester 

& Taylor 1977) B12 = 3.2 x WVpP, we then have 

i?wind = llM^^Rl ohms. (13) 

For the potential drops of inner gaps (ICG and lAG), when the spark provides necessary 
charges to close the pulsar circuit, they should correspond to resistors (i?iAG and -Ricd shown 
in Fig.2). When there is no sparks, the gap grows and they can be described by capacitors 
(CiAG and CicG, shown in Fig.2). 

Let's analyze the circuit in Fig.2. Although the electric power has fixed potential supply, the 
current is changing due to the inner gap sparks. In this sense, the equivalent circuit description 
in this paper is not simply a DC circuit model. Because of the erratic sparking, the resistance 
i?RS could be as the sum of many sinusoidcil functions of time, R^,s — ^^^o-^'^^ 

sin riLot. The 

electric current between arbitrary points M and N in circuit could also be in this form, 7mn = 

I]5^q/„ sin nujt. 

According to the Kirchhoff 's current and voltage laws, the complex impedances of parallel 
connection circuits composed of i?icG and CicG) Riag and Ciagj -Rwind and Caow are zicG = 

RiCg/ (1 + i^RlCcClCG) ■ ZlAG = Riag/ (1 + it^-RlAGClAc). ^wind = -Rmdnd/ (1 + iW-RwindCflomO + 

iw-Lflow) respectively, where i = a/— T, oj the angular frequency of electric current modulation. 
Defining z' = zicg + ^iag + -^wind, one obtains the total complex impedance to be 

z' 

Ztotal = rr—. — 77; ohms. (14) 

Let u) = cjo when |a;totai| is the smallest value. In this case, the potential drop between inner 
vacuum gap is the highest. It is possible that there exists an oscillation with time scale 
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in the circuit. We expect that some of the variations of radio intensity in different timescales 
could be hints of such circuit osciUations. 

In case of Cjcg = Qag = ^flow = C'flow = 0, one has the total impedance z = i?icG + 
-RiAG + -Rwind- The physical meaning of this result is of the DC circuit model (Shibata 1991), 
where i?rot — 27rrpCpGJ ■ K;ap7 which results in the deduction of pulsar magnetic fields B ~ 

^c?IPP/{-K^R^) (e.g., Manchester & Taylor 1977), is usually assumed (e.g., Xu & Qiao 2001), 
with the Goldreich- Julian density poj. 



3 PULSARS CHARGED ELECTRICALLY? 

If charged particles distribute as the Goldreich- Julian density, pgj i they will be in balance about 
electrostatic force. The equivalent "Possion" equation in comoving frame then is (e.g., Beskin, 
Gurevich & Istomin 1993) 

V -E^A^ip- PGi). (15) 

As mentioned by RS75, the electric field, £^gJ: on star surface due to the lack of charge density 
respective to the Goldreich- Julian density (e.g., for vacuum outside the star, p = 0) is normal 
to star surface. The solution of RS75 (see its appendix 1.6) for the electric field on the stellar 
surface is Eg = ~2VLBh/c < 0, which is equivalent to that of choosing the potential of field line 
a to be the same one of the ISM, 0^ = 0. From equation H15|l . one can also find generally the 
surface electric field Eqj < in both Regions / and //, as shown in Fig. 3. 

When the potential of the line b is zero, as mentioned above, the vector direction of electric 
field El within Region lis inward but that in Region 11, En, is outward. How can one understand 
consistently this picture? Why should this be reasonable if one choose 0b = 0, rather than 
(pA = 0? The answer could be that there must be positive charges on the star surface which 
increases the potential of star. These charges provide a monopole electric field Ej^o, which 
combined with Eqj to form the total electric field, as shown Ei and -En in Fig. 3. 

The electricity induced by the charges on the star surface is so high that the electric field is 
reversed across Region //when electric field of star self and the field caused by charged particles 
are combined. Therefore, the current can flow out and come back in magnetosphere to close the 
pulsar's generator circuit. Reversely, the increased electric field supports our assumption that 
the potential of line b is zero. 

Provided that two conditions are satisfied (1, the star has positive charges on the surface; 
2, the potential of the line b is zero), the potential increase of star relative to that of ISM could 
then be estimated to be order of the potential drop between A and B (from equation 2) 

T/ W> W> 3><10l6i?|Bi2, . 2^ . 2. ^ 3 X 10^^RlBi2 ,^ , . 

^star 0A - 0B = 7z (sm - sm 9b) ~ 7^ volts. (16) 

B B'^ 

From equations Q and 116|) . the total charges on stellar surface could then be 
„ „ 3 X 10 ^R\B\i 

Q = Kitar X Cstar ~ COulombs. (17) 

The collapse of an evolved massive star should form temporally a rotating magnetized compact 
star, and then become a black hole. The black hole could be charged too, with a quantity to 

be proportional to B/P^. As magnetic field B ^ \/ Pp, one has then Q ~ Jp/P^ ^ e]/'^ 



-'lot ' 

where the rotation energy loss rate E^ot Observationally, the X-ray luminosity could be 

a function of the spin-down energy loss for all rotation-powered pulsars, Lx ~ E^ot (Becker & 
Triimper 1997); but the 7-ray luminosity Lj ^ E^.^^ (Thompson 2003). One can also address 
that energetic 7-ray luminosity is proportional to the electric charge Q, rather than to /?iot- 
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Fig. 3 An illustration of the total electric field E as a. function of polar angle 0. Eqj 
is due to the charge-separation in pulsar magnetosphere, provided that the potential 
of line a is the same of the ISM. E^no is induced by the charges of star. The total 
electric field (middle solid thick curve), Ei and En within the corresponding regions, 
could be the combination of Eqj (lower solid thin curve) and E^o (upper solid thin 
straight line). 



One thinks conventionally, according to equation H15|) . that no acceleration (e.g., E'n — 0) 
occurs if p = pGJ in pulsar magnetosphere. We note, however, that this conclusion is valid only 
if no solenoidal force field appears. In other words, i^y ^ though p = pcj if one adds any 
solenoidal force field in the magnetosphere. A charged pulsar contributes a solenoidal electric 
field, which results in an acceleration near the pulsar (to damp as 1/r^). In the close field line 
region, this electric field causes a re-distribution so that E^^ = 0. In the open field line region, 
this field accelerates particles (for parallel pulsars, to accelerate negative particles in Region 
I, but positive particles in Region II; see Fig.l). Certainly, extra acceleration due to p ^ poj 
exists too. As demonstrated in Fig. 3, a very large acceleration may exist near the last open 
field lines, which could be favorable for the high energy emission in the caustic model (Dyks 
& Rudak 2003). An outer gap may not be possible if particles can flow out freely either from 
the surface (for negligible binding energy) or from the pair-formation-front (for enough binding 
energy) of a charged pulsar. 



Pulsars could be bare strange stars, some of them could be of low-mass (Xu 2005). Due to 
the color self-confinement of quark matter, the density of low-mass bare quark star is roughly 
homogeneous, and its mass would be 



where the bag constant B = 60i?60 McV fm ^ (i.e. 1.07 x lO^^'g cm •^). For a star with pure 
dipole magnetic field and a uniformly magnetized sphere, the magnetic moment is 



4 LOW-MASS BARE STRANGE STARS 




(18) 




(19) 
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If the magnetized momentum per unit mass is a constant fj,^ — (10 ~ 10 ^)G • cm^ • g ^, the 
magnetic moment then is (Xu 2005) 

fi = n^M. (20) 
Combing equations H18() - (|20|l . one can obtain the magnetic field stength 

B = 1.8 X lO^^VmSeoMo. (21) 

Therefore, for a low-mass bare strange star, if the values of its period P, radius R (or the 
mass Mqg), and the polar magnetic field B (or the parameter /im) are all known, the total 
charges on the surface could be obtained by equation (|17|l . 



4.1 Evidence for low- mass millisecond pulsars? 

Pulsar's radius is assumed as a constant in above sections. In fact, the radius should be a 
variable in different models (of either normal neutron star or strange quark star). The radius 
of bare strange stars could be as small as a few kilometers (even a few meters) . Could one find 
any observational hints about the star radius? It is suggested that normal pulsars could be bare 
strange stars with solar masses, whereas millisecond pulsars are of low masses (Xu 2005). Can 
we show evidence for low-mass millisecond pulsars in the bare strange star model for pulsars? 
These are investigated, based on the observational pulsar data^. There are 1126 pulsars with 
known P, P and radio luminosity L1400 (mJy kpc^) at 1400 MHz simultaneously. The numbers 
of millisecond and normal radio pulsars are 35 and 1091, respectively, if the dividing line of 
them is P = 15 ms. Here we assign the 35 millisecond radio pulsars (P < 15 ms) as Sample I 
and the 1091 normal radio pulsars (P > 15 ms) as Sample II. 

Assuming the pulsar radius is a variable, from equation H17(l . one can obtain pulsar's charges 

Q ^ R'^^ ^ R'^i^Y^^ since the magnetic field strength B ~ \/ PP. Defining ( = one 
comes to 

logQ~41ogP+ilogC. (22) 

At the same time, the rotation energy loss rate E'lot — I^^ ^ R^C: where rotational inertia of 
star / ~ ^ME? ~ (the star is assumed to be a homogeneous rigid sphere), i.e., 

log4ot-^51ogP + logC. (23) 

The correlation between L1400 and Q and the normalized statistical distribution of L1400/C 
are shown in Fig. 4. For showing the evidence for low-mass millisecond pulsars by this relation 
and distribution, we firstly give the best fit line by the least square method for larger Sample 
II. Secondly, we assign the slope of fitting line for the smaller Sample I to be the same as that 
of Sample II. Thirdly, we study the relation of L ^ R with two assumptions, respectively, that 
are the luminosity L = L{Q) (i.e., L is only a function of Q) and L = L(Ej-at) (i-e., L is only a 
function of -Ejot)- Finally, through comparing the intercepts of two fit lines, we can obtain the 
ratio of radius between Sample I and Sample II. 

From left panel in Fig. 4, we find that the best fit line for Sample II are 

logLi4oo = 2.86 + 0.11 log C (24) 

Assigning the same slope for Sample I, we find that the intercept of fit line for this sample is 
1.56. 



* http://www.atnf.csiro.au/research/pulsar/psrcat/ 
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Fig. 4 Left: Observed luminosity L1400 (mJy kpc^) at 1400 MHz as a function of 
C for 1126 radio pulsars, where ( = The solid squares represent the millisecond 
radio pulsars (Sample T) with P < 15 ms while the empty circles are the normal ones 
{Sample IT) with period P > 15 ms. The solid line is the best fit result with the least 
square method for Sample II, and the dash line is the fit for Sample / assigned the slope 
to be the same as that of Sample II. Right: The normalized statistical distribution 
of ii4oo/C for Sample I (dash step line) and Sample II (solid step line). The vale of 
-^^i40o/C is only a function of pulsar radius according to Ea. l|22l24|l . 



If L ~ Qj from equations (22) and (24), the ratio of radii Rn for Sample II and Ri for 
Sample I at 1400 MHz then should be 41og^ = 2.86 - 1.56 = 1.30. Then one can obtain 

Rii « 2. IIP/. In the same way, if L ^ Prot, from equations (23) and (24), we can obtain 
Rii w 1.82P/. From right panel in Fig. 4, it is evident that the values of P1400/C; a function 
of pulsar radius, distribute with two peaks; the higher one (i.e., larger radius) for normal pulsars 
whereas the lower one (i.e., smaller radius) for millisecond pulsars. 

Summarily, the radii of Sample II are always larger than those of Sample I, which gives an 
evidence that the radii of millisecond pulsars are smaller and the masses are lower than those of 
normal radio pulsars. Especially, from the right panel in Fig. 4, we can find that the normalized 
distributions of P1400/C for Sample / and Sample //are almost clustered at about 13 and 15. 
This two-peak structure gives a crude same result and thereby supports again the evidence for 
low-mass millisecond pulsars. If one thinks that normal pulsars are bare strange stars with mass 
~ Mq, millisecond pulsars could be bare strange stars with mass ~ (1/2)'^ Mq. 



5 CONCLUSIONS 

Assuming that the magnetosphere of pulsar has a global current which starts from the star, 
runs through the inner core gap, the wind, the outer gap and the inner annular gap, and returns 
to the star, we study the circuit characteristics of four elements: the electromotive source, the 
inner core gap, the inner annular gap, and the outer gap and wind. It is emphasized that the 
potential of critical field lines equals to that of interstellar medium. We find, in this case, that 
the pulsar whose rotation axis and magnetic dipole axis are parallel should be positively charged 
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for the pulsar's generator circuit to be closed. The current flows out through the light cylinder 
and then flow to the stellar surface along the open magnetic field line. 

There are five independent parameters which can describe completely the dynamics of a 
pulsar magnetosphere assuming a dipole field and a uniform density of the star. They are the 
radius i?, the mass M(~ WirBR^ /3 in case of bare strange stars, with B the bag constant), 
the magnetic strength B (or magnetic moment /i ^ BR^), the period P, and inclination angle 
a. Typical parameters for radio pulsars are: R ^ 10^ cm, M ~ IAMq, P ~ (10^'^ — 1) s, 
B ~ 10^"^^ G. No solid observational evidence shows these parameters are really typical. In 
case that pulsars are bare strange stars (probably with low masses), some of the parameters 
above may not be representative, and the independent parameters could be related each other 
(Xu 2005). The statistics between the radio luminosity and pulsar's electric charge (or the 
spindown power) may hint that millisecond pulsars could be low-mass bare strange stars (with 
masses of a few O.IMq). 
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